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Abstract. - I introduce an ansatz for the exclusion statistics parameters of fractional exclu- 
sion statistics (FES) systems and I apply it to calculate the statistical distribution of particles 
from both, bosonic and fermionic perspectives. Then, to check the applicability of the ansatz, I 
calculate the FES parameters in three well-known models: in a Fermi liquid type of system, a 
one-dimensional quantum systems described in the thermodynamic Bethe ansatz and quasiparticle 
excitations in the fractional quantum Hall (FQH) systems. The FES parameters of the first two 
models satisfy the ansatz, whereas those of the third model, although close to the form given by 
the ansatz, represent an exception. With this ocasion I also show that the general properties of the 
FES parameters, deduced elsewhere (EPL 87, 60009, 2009), are satisfied also by the parameters 
of the FQH liquid. 
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Introduction. — Haldane's concept of fractional ex- 
clusion statistics (FES) [1] have been recently amended 
in a series of publications [2-5]. In these publications I 
showed that in the original formulation of FES some ba- 
sic properties have been overlooked. The situation was 
corrected first by introducing a conjecture [2] and then 
by deducing the general, basic properties of the FES pa- 
rameters [5]; it turned out that the conjecture of Ref. [2] 
is just a special case of the general conditions deduced in 
Ref. [5] , which allows one to write down an explicitely con- 
sistent system of equations for the statistical distribution 
of particles in a FES system. 

Now, that the general properties of the FES parameters 
are deduced [5], the conjecture of Ref. [2] looses its status 
and becomes simply an ansatz. This ansatz seems to be 
quite general and applies to most of the macroscopic (i.e. 
quasicontinuous) systems. Nevertheless, by the end of this 
letter I will show an exception. 

In this letter, using the general properties of the di- 
rect exclusion statistics parameters (see eq. [2c] below) I 
propose an even more restrictive (and therefore more con- 
veninet for applications) form of this ansatz, by specifying 
also the form of the direct FES parameters. This new 
ansatz allows me to write the system of equations for the 
statistical distribution of particles in a more clear form 



and to single-out the direct exclusion statistics parame- 
ters which are most often used in FES calculations. 

In general, the FES equations for the particle distribu- 
tion are used in the bosonic formulation. Here I will write 
and use these equations in both, bosonic and fermionic 
formulations. This allows for direct application of the for- 
malism in systems of cither bosons or fermions. 

In the end I will give three examples of systems in which 
FES its manifesting. After I calculate their parameters, 
I show that two of them satisfy the ansatz, whereas the 
third, although quite similar, constitues an exception. 

The general properties of the exclusion statistics 
parameters. — Let us assume that we have a system 
of particles which we divide into the species indexed by 
i = 0, 1, . . ., each of the species containing Ni particles 
and having Gi "available single-particle states" . Then the 
number of microscopic configurations for species i is 



A,!G,! 



(1) 



To introduce FES into the system, we define the exclusion 
statistics parameters, ctij , such that at the variation of Nj 
by SNj, the number Gi changes by SGi = —ctijSNj. The 
diagonal elements, an, are called direct exclusion statistics 
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parameters, whereas the nondiagonal ones, , i ^ j, are 
called mutual exclusion statistics parameters. 

In Ref. [5] I showed that if one of the particle species, 
say species j, is divided into the subspecies jo,ji,--., then 
the new statistics coefficients should satisfy the relations 



a 



'^3031 



for any i, i ^ j (2a) 
for any i, i ^ j (2b) 

. = ... (2c) 



Here I am using the notations ctij, like in Refs. [2-5], to 
make the difference between the "extensive" and "inten- 
sive" FES parameters that are going to be defined below. 

The conjecture introduced in Ref. [2] stated that in a 
macroscopic physical system described by FES there is a 
division of the system into species, say {{Gi, Ni)}i=o^i^,,,, 
so that no matter how we divide further these species 
into subspecies, the mutual exclusion statistics parame- 
ters are always proportional to the dimension of the space 
on which they act. Concretely, this means that for any 
i, j, with i j, we can write 5y = QijGi, where are 
constants that depend on the species i and j, and at any 
further division, say species i is divided into the subspecies 
io,ii, ■ ■ •, the new mutual exclusion statistics parameters 
satisfy cti^j = aijGi^^ for any i ^ j [2]. The parameters 
oiij are then called the extensive parameters and the pa- 
rameters aij are called the intensive parameters. It is easy 
to check that the extensive parameters satisfy the general 
conditions (Pa|) and (j2bp . 

In this paper I will extend the ansatz to the direct FES 
paramteres, which should satisfy (|^ . 

Ansatz for the FES parameters. — We can make 
the form of the exclusion statistics parameters that satisfy 
Eqs. ^ more specific and easier to apply to FES calcula- 
tions if we decompose fiy into a sum of two different types 
of parameters, dtf^ and dif\ by the relation 



statistics in the system. In such a case the thermodynamic 
calculations simplify considerably. Note also that the ideal 

~ (s) 

Fermi gas corresponds to = 1 for any i. 

Particle population in the bosonic formulation. — 

Let us now deduce the equations for the particle popula- 
tion. To avoid unphysical (negative or divergent weights) 
I write the number of microscopic configurations as [2] 

„ (G, + iV, - 1 + (1 - a^;^)5N, - E, a^j5N,)\ 



3 



n 



{N, + 5N,)\{G^ - 1 - af>5N, - ^ 
[G, + iV, + (1 - a^'^)5N, - Q a^5Nj\\ 
[N, + 5N{)\{G^ - a^'^SN, ~ G,; E, a,,5Nj)\ 



which then I plug into the expression for the grandcanon- 
ical partition hmction. 



^ I¥({A^,})exp 



(7) 



where /3 = l/fc^T, T is the temperature of the system, 
whereas /^i and ti are the chemical potential and the 
single-particle energy for the particles of species i. Max- 
imizing Z with respect to the populations ni = Ni/Gi, I 
obtain the system of equations 



PilJ-i 



In- 



(8) 



Notice that Eq. jS]) is similar to Eq. (18) of Ref. [2], 
only that by singling out the coefficients ai we could 
extend the summation on the r.h.s. to include also the 
terms i = j. This makes the application of Eq. ([8]) more 
straightforward than the one in Ref. [2] . 

If aij = for any i and j, we recover the typical formulas 
for the calculation of particle population without mutual 
exclusion statistics [6] , 



a. 



(3) 



a) 



(9) 



The parameters a\f 



are the "extensive" 



above, only that now, by separating a, 
the condition 

,G 



~(e) _ 



ones discussed 
we can extend 

(4) 



-is) 



also to the case i = j. 

The additional parameters, 
one species of particles and are not extensive. If we split 
the species i into the sub-species ioj^ii ■ • -j then by Eqs. 



always refer to only 



with w and C defined by 

nf^' (C)[l + ^i'(C)]'-"-''' = C ^ exp[/3(e, - ^i,)] (10) 

In the quasicontinuous limit, in a phase-space spanned 
by the single-particle states of quantum numbers k (k is 
not necessary the wave-number), of density of states cr(k), 
Eq. ([5]) transforms into 



cl, ([3]) and the extcnsivity property of otfi , we obtain /3(/Xk — ek)+lii 



[l+nkP^V 



(5) 



Typically, in the literature we find exclusion statistics 
parameters of the (s) type (see e.g. [6-13]). Therefore 
in general onj = for any i ^ j, so there is no mutual 



= / (T(k')ln[l+nk']ak'kdk'. 

(11) 

Equation (|lip is similar to eq. (19) of Ref. [2], if we iden- 
tify, say 01^^ with cikk of [2]. The identification is natural, 
since as the interval Gi becomes smaller, an converges to 



a\^' , which stays constant, whereas oil. 



decreases to zero. 



p-2 



An ansatz for the exclusion statistics parameters in macroscopic physical systems described by fractional exclusion statistics! 



Particle population in the fermionic formulation. 
— Formula ([T]) represents the number of configurations of 
particles of species i in the bosonic formulation of FES [1] . 
This description is not the most convenient for example 
when one describes systems of (interacting) fermions in 
the FES formalism, since in such a case Gi represents 
the difference between the number of single-particle states 
and the number of fermions. Therefore in such cases it 
is easier to work directly with the total number of states, 
Ti = Gi + Ni — 1, in the fermionic description. Although 
the two descriptions are equivalent, let me write down the 
system of equations for the n^s in the fermionic descrip- 
tion. For this, I write first the number of configurations 
for the species i. 



In Eqs. (Uni) and pT)) (i = 0, 1, . . .) are single-particle 
energies, rii is the population of the state i, and Vij repre- 
sent the interaction energy between a particle on the state 
i and a particle on the state j. 

Going to the quasi-continuous limit, assuming that the 
single particle energy spectrum has the density of states 
a{e) and that the interaction energy depends only on the 
energies, we replace the indices i and j by the energies e 
and e' to write 



(/) _ 



T,} 



(12) 



We add again small perturbations to the particle numbers 
and write 



1 



{N, + 5N,)\ 



I plug (fT3|) into the expression ^ for Z and, by maxi- 
mization, 1 get the equations for the particle population, 
fi = N,/T,, 

Introducing the density of states a'^^^ (k) I write Eq. 
in the quasicontinuous limit. 



^(A*k - Ek) + In 



/k 



= - / aW)(k')ln[l-/k']ak'kdk'. 



(15) 



Applications. — Let me now analyse three interact- 
ing particle system models in which FES is manifesting 
and compare their FES parameters with the ansatz pro- 
posed here. 

FES in a system with Fermi liquid type of interaction. 
I take again the model of Ref. [3] , which is a generalization 
of the Murthy and Shankar model [7] , widely used in FES 
[7-10, 12, 14, 15]. In this model the total energy of the 
system, 

i ij 

is splitted into the quasiparticle energies as E 
with 

i-l 



1 



(16) 



(17) 



e + 



V{e,e')a{e')n{€')de'. 



(18) 



(13) 



In what follows I shall assume that the function e(e) is 
bijective and therefore I shall use interchangebly, whenever 
necessary, both e(e) and e(e). 

The FES is manifested in this case in the quasiparticle 
energies, e. In order to describe FES and to calculate its 
parameters, we split the e axis into small intervals. In 
general we shal denote such an interval by Se and by this 
notation we shall assume that it contains the quasiparticle 
energy level e. 

From Ref. [3] we can directly identify the direct exclu- 
sion statistics parameters of (s) type. 



a): 



Vie,e)a{e) 



1 + 



dVie,e') 



a{e')n{e')de' 



(19) 



where e = e(e). 

The mutual exclusion statistics parameters are [3] 



aseSe 



V{e,e,) + fCe,e^) 
''l + i;^^'j{e'He')de' 



daje) 
de 



(20) 



where e = e(e), = e(ei), e G Se, and q £ 6ii. For ([20|) . 
the quasiparticles are inserted into the interval Sii while 
the variation of the number of states is observed in the 
interval Se. The function f{e,ei) is 



dV{e,e') 



a(6')n(e') 



Se' 



de', (21) 



where by p{e) = cr(e)n(e) I denoted the particle density 
along the e axis and the notation [i5eV'^P(£i)]{p(g)} repre- 
sents the functional derivative of e' with respect to the 
particle density at energy e^, when we keep fix e(e'). Since 
the number of states in the interval is a{e)5e, the coef- 
ficients agg, are 



F(6, £,)+/(£,£.) 



de 



a{e')n{e') de' 



dlogcr(e) 



de 



(22) 



(e) 

Now we are left with the calculation of a~i . For this 
we first note that Eqs. ([20|) and ([22|) are valid for any two 
disjoint intervals, so let us divide the interval Se into the 
subintervals Sii, i = 0,1,..., M, of dimensions a{ei)Sei, 
where we always maintain the convention that e,; belongs 
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to the interval Sii. We pick the energy level Cfc from the 
interval ^e^ and apply Eq. (j2cp : 



asiSe = o-(eo)(5eoajoj^ +. . . + asi^sik + - ■ ■ + <^{iM)5eMaiMek 

(23) 

Making the interval ^e^ small enough as compared to 5e 
and assuming that 5i is also small, so that we can use 
the approximations agj.g; = agg and (T(efc) = for any 
k,l = 1, . . . , Af, we can simplify Eq. ([23]) to write the 
general expression 



~{s) _ 



aud[e)5i+5{l-h)a'^^\ (24) 



which has the form of the ansatz ([3|). For Eq. (f24|) . note 
that /{e,e) = 0. 

In the simplified models of Refs. [7-10, 12, 14, 15], only 
the direct exclusion statistics parameters have been used 
and agg; was identically zero for any e and e^. We observe 
now, from eq. (j20p . that this happens whenever da/de — 
0. 

Having all the exclusion statistics parameters calcu- 
lated, one can in principle apply Eq. pT|) or ([TS]), de- 
pending on the type of particles we have in the system, to 
calculate the particle distribution. 

FES in ID quantum gas in the thermodynamic Bethe 
ansatz. The ID gas of quantum particles in the thermo- 
dynamic Bethe ansatz (TBA) have been analysed before 
(see e.g. [9, 16, 17]) and is recognized in general as being a 
system which can be described by FES. The only reason 
for which I discuss it again here is to show that its FES 
parameters are indeed of the type ([3]) and also because in 
general a confusion is made in the literature between the 
intensive a^- and extensive 5^ parameters and this has to 
be clarified. 

Therefore let's consider the typical gas of N spinless 
particles, bosons or fermions, on a ring of circumference 
L. We assume that the system is nondiffractive [18] and 
the asymptotic particle wavenumber, k, is determined by 
the equation [19] 



Lk^J29{k^k')^Lk^°\ 



(25) 



where fc*-"-* = 27r/(fc)/L is the free-particle wavenumber, 
I{k) an integer that depends on k, and d{k — k') is the 
phase-shift due to the interaction. 

To simplify the notations and to be also in accordance 
with Refs. [18, 19] we set the units so that h — m = 1, 
where m is the mass of the particle. In these units the 
total number of particles, momentum and energy of the 
system are 



TV-^l, P = Y.k, and£;=i^/c2 



(26) 



For large systems we transform the summations into in- 
tegrals and define the densities of states, cr(fc) and cro(fc(°^), 
by the relations 

D{Sk) EE cr{k)Sk and D{Sk^°')) = (To(fc("^)(5fc(°) (27) 



where L»((5fc) and ^'(JfcW) are the numbers of states in the 
small intervals 6k and respectively. If Sk and Sk^^^ 

are related by Eq. then D{Sk) = D{Sk^°^). Obvi- 

ously, iTo(fc^*''') = L/{2'k) (if we impose periodic boundary 
conditions on k^°^) [18,19]. 

The populations of the single particle levels, no{k'^*^^) 
and n(fc), are defined as 

N{6k) = n{k)a{k)6k = N{Sk'-°^) = noik'^°'^)(7o{k^°^)Sk'-°^ 

(28) 

where N{dk) = N{6k^^^) is the number of particles in the 
interval 5k or Jfc^''^ with 5k and related by Eq. psj) . 
Moreover, since both, the number of particles and the 
number energy levels, are the same in the 5k and 
intervals, we have the identity n{k) = no[kS^^{k)]. 

In accordance with the notations in the literature 
[16, 19] I introduce here also the particle density, p{k) = 
a{k)n{k)/L = N{5k)l{L5k). 

In the new notations, Eq. ([^5]) becomes a self-consistent 
equation for k, 



k^k'^°\k)+ / e{k-k')p{k')dk', 



(29) 



from which we can calculate [19] 



1 - / e'{k - k')p{k')dk', (30) 



dk 



where 9'{k) = de{k)/dk. Plugging Eq. ^ into (E?]), we 
get the density of states 

a{k) = ^1^1- J e'{k^k')p{k')dk'Y (31) 

The FES is manifesting in the system because of the 
dependence of a on p: the variation of p{ki) (changing 
the number of particles at the level ki) produces, in prin- 
ciple, a change of the density of states cr(fc), at any k. 
To determine the coefficients of the exclusion statistics we 
calculate the variation of cr(fc) at the variation of p{ki), 
i.e. we calculate the functional derivative 



L5p{k,) 2tt ^ " 



(32) 



Therefore if we split now the k axis into small intervals, 
with 5k the interval around k and 5ki the interval around 
ki, then the variation of N{5ki) by 5N{6ki) produces a 
variation of D{5k) equal to 

5D{5k) ^ /r^,t\ (5fc = -^0'{k - ki)5k5N(6k,) (33) 
Ldp(ki) ZTT 

From the FES formula, 5D{5k) — —askSk^5N{5ki), we 
obtain imediately 



askSk, = -^d'{k - ki)5k 

ZTT 



(34a) 
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and therefore the exclusion statistics parameter askSki is 
proportional to the dimension of the space on which it acts, 
D{Sk). 

Comparing Eq. (|34ap with Eq. ^ we get 



a-kki 



9'{k - kj) 
2TTa{k) 



(34b) 



From Eqs. ([M]) all the TBA thermodynamics follows in 
general, by direct application of the formalism presented 
before, in either bosonic or fermionic formulations, as we 
shall see imediately. 

If the particles are bosons, we plug auki into setting 
di^^\ we obtain after some simplifications 

keT 



e{k) = ef-^,- 



27r 



\og[l~e-'^'^^"^]e'{k'-k)dk'. (35) 



In Eq. (I35|) I used the notation e^"'' = k'^ /2 and I defined 
the quasiparticle energy, e(fc), by = {exp[/3e(fc)] — 

If the particles are fermions, we define the quasiparticle 
energy by exp[— /3e(/c)] = nfe/(l — Uk) and we plug Eq. 
pib]) into Eq. dH]). In this way we recover imediately the 
TBA equation. 



e(fc) 



.(0) 



{k) 



keT 
27r 



log[l + e-'^'=('^')]6l'(fc'-fc) dk', 

(36) 

similarly to Eq. ((35)) . 

For the delta function interaction potential between the 
particles, V{x) = 2c5[x) {x being the distance between 
the particles), the phase shift is 9{k) = — 2 arctan(fc/c), 
which gives 



a-kki 



and therefore 



-(e) 
'^SkSk, 



na{k) c? + {k ~ kiY 



5k 



TT c2 + (fc - fci)2 ' 



while 



(37) 



(38) 



(39) 



for any k. 

If the particle-particle interaction is V{x) = A(A — l)/x^, 
then the phase shift is 6{k) = 7r(A — l)sgn(fc), where sgn 
is the sign function. From this we obtain 



O-kki 



(A-1) 
<j{k) 



S{k-h). 



(40) 



Note that in this case we do not have any "extensive" 
mutual exclusion statistics parameters and therefore we 
may write askSk, = ct-^^'^SkSk, = (A - l)5skSki- 
In both cases we recover the ansatz Q. 



Fractional Quantum Hall Effect. Another system 
which is traditionally related to the FES is the fractional 
quantum Hall effect (FQHE) [1,6,20,21]. 

In a Laughlin 1/m-liquid, with m an odd integer, at 
any finite temperature there are quasiparticle vortex-like 
excitations, Nj^ and corresponding to quasi electrons 
and quasi holes. The numbers of quasi excitations are 
related to the number of flux quanta in the system, = 
eBA/hc (where B is the magnetic flux and A is the area 
of the sample), and electron number, N^, by the relation 
[20,21], 



(41) 



For single particle ocupancy (only one quasi-excitation in 
the system), the number of available states for each of 
these types of excitations is G_ = G+ = N^, while for 
general iV+ and we have 



1 



(42a) 



G_ = —Na, - a^+N+ - a^^N^ (42b) 
m 

where {i,j — +, — ) are the FES parameters. Although 
maybe there is still no consensus regarding the values of 
aij, which differ for the different liquid models used to 
describe the FQHE (see e.g. [21]), for concreteness I shall 
adopt here the bosonic vortex scheme [20-22] , with 

a++ = 2 — 1/m, = l/m, (43a) 

a |_ = —l/m, a = l/m, (43b) 

alhough this is not important for our discussion. 

The point is that although the fractional quantum Hall 
liquid (HQHL) is a macroscopic system, apparently the 
FES parameters of this system do not obey the general 
relations ([2]): there are only two species of quasiparticles- 
the quasi-electrons and the quasi-holes, with degenerate 
energy levels-and the a's are fixed by (|43)) . so we cannot 
(aparently) split the species into further subspecies. Still, 
relations ([2]) are deduced on very general grounds, so they 
should be valid. 

We have a puzzle. 

The solution of this puzzle is straightforward and may 
be obtained only by macroscopic considerations. For this 
we observe that FQHL being a macroscopic system, one 
can always divide its area. A, into smaller areas, Ai, 
i = 0, 1, . . ., and in each of the smaller subsystems the 
FQHL has the same properties, but with scalled number 
of electrons, Ne-, flux quanta, N^., quasi-electrons, N+., 
quasi-holes, N^., and G±. ~ G± ■ Ai/A. In equilibrium, 
iVe, = N, ■ A,/ A, N^^ ^' ■ AJA, iV+. = N+ ■ AJA 
and N^. — ■ Ai/A. In each of these subsystems, the 
same relations, (j42|) and ([43]) . are valid for the quanti- 
ties Ne,, N^,, N+^, N^., and G±, = G± • Ai/A. Notice 
with this occasion that there is no "mutual" statistics be- 
tween different sub-systems, i and j [i ^ j). Therefore, 
by continuing to split A into samller and smaller areas, 
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we may eventually end-up with a coarse-graied surface of 
elemetary areas, SA{x,y), where the x and y are the two- 
dimensional coordinates on the surface of the FQHL, and 
with the FES parameters, aij[SA(x,y)] ~ — ) that 
act always on the same elementary area, {dA{x,y)). 

The puzzle is only aparent. Since the system is macro- 
scopic, therefore extensive, the quantities ay are local and 
rigorously should be written as (r, r') = aijS{r — r'), for 
any i, j = +, — , where by r and r' I denoted two position 
vectors {x,y) and {x',y') in the plane. This form of FES 
parameters obey the genreal relations ([2]), as they should, 
but constitue an exception to the ansatz ([3]) proposed in 
the beggining, due to the fact that the off-diagonal ele- 
ments a_| and a |_ are not extensive, but are also pro- 
portional to (5(r — r'). 

In the reverse process, if we "glue" together all the el- 
ementary areas, SA{r), we reobtain the original system, 
of area A, two species of particles, iV+ and N^, and the 
overall FES paramters (I43p . 

Conclusions. — In this paper I introduced an ansatz 
for the FES parameters and with it I calculated the sta- 
tistical distribution of particles from both, bosonic and 
fermionic perspectives. This ansatz allowed me to write 
the system of equations for the statistical distribution of 
particles in a more clear form and to single-out the direct 
exclusion statistics parameters which are most often used 
in FES calculations 

Then I took three examples: a Fermi liquid type of sys- 
tem, a one-dimensional integrable quantum system and a 
fractional quantum Hall (FQH) system. I calculated the 
FES parameters for these systems and I showed that those 
of the first two systems obey the ansatz proposed here, 
whereas those of the third one do not. With this ocasion 
I also showed that if one takes properly in to account the 
extensivity of a FQH system, then its FES paramters also 
obey the general conditions deduced in Ref. [5]. 
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